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Abstract 



We propose a new, twistor string theory inspired formalism to calculate loop ampli- 
tudes in M = 4 super Yang-Mills theory. In this approach, maximal helicity violating 
(MHV) tree amplitudes of M = 4 super Yang-Mills are used as vertices, using an off- 
shell prescription introduced by Cachazo, Svrcek and Witten, and combined into effective 
diagrams that incorporate large numbers of conventional Feynman diagrams. As an ex- 
ample, we apply this formalism to the particular class of supersymmetric MHV one-loop 
scattering amplitudes with an arbitrary number of external legs in M = 4 super Yang- 
Mills. Remarkably, our approach naturally leads to a representation of the amplitudes 
as dispersion integrals, which we evaluate exactly. This yields a new, simplified form for 
the MHV amplitudes, which is equivalent to the expressions obtained previously by Bern, 
Dixon, Dunbar and Kosower using the cut-constructibility approach. 
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1 Introduction 



It has been proposed recently [1] that J\f = 4 super Yang-Mills (SYM) has a formulation 
as a topological string theory in super twistor space. This proposal, and related formula- 
tions, have since been actively investigated (see [l]-[25]). The original maximal helicity 
violating (MHV) tree amplitudes given in [26, 27] take a surprisingly simple form, and 
the tree amplitudes in the gauge theory have now been understood in the context of this 
topological string theory. Inspired by the twistor approach, the MHV amplitudes were 
used as vertices by Cachazo, Svrcek and Witten (CSW) in [6], together with a particular 
off-shell prescription, in order to present a new and extremely efficient method of deriving 
more general, non-MHV tree amplitudes. This has been further explored in [7,13,20-24]. 
The CSW method of constructing gauge theory tree amplitudes corresponds in the twistor 
picture to completely disconnected instantons of degree one, linked by twistor space prop- 
agators, whereas the original approach of [1] used connected higher degree instantons. The 
picture emerging is that there appears to be a number of alternative ways to compute 
tree amplitudes, which are related by the degeneration of higher degree curves into curves 
of lower degree [2,4, 10, 14,21]. 

The study of quantum amplitudes in Af = 4 SYM has led to many interesting results 
- see [28]- [38] and references therein. Much of this progress was achieved thanks to 
powerful unitarity-based methods [30,31,33,34,38], whereby properties of loop amplitudes 
are derived from cutting rules and tree amplitudes. One of the main virtues of this 
approach lies in the fact that the tree amplitudes on both sides of the cut are evaluated 
on-shell, and can moreover be simplified before constructing the loop amplitude. This 
is referred to as the "cut-constructibility" approach. The simplest application concerns 
the calculation of one-loop MHV diagrams (see [32,33] for reviews). Here one can derive 
explicit expressions for the amplitudes at one loop using the known results for the tree- 
level MHV amplitudes [29]. Recent work [36,37] has also uncovered intriguing recursion 
relations which link two-loop amplitudes to those at one loop, and led to speculation about 
the existence of more general cross-order relations between L — 1 and L loop amplitudes. 

So far, twistor-inspired techniques have successfully reproduced known MHV and non- 
MHV tree amplitudes and inspired efficient new techniques to calculate new tree ampli- 
tudes. One expects that this success will carry over to perturbative calculations at the 
loop level. Whilst some initial comments on loop diagrams were already made in [1] and, 
more recently, the twistor space structure of one-loop diagrams in M = 4 SYM has been 
analysed in [25], it has so far not been clear how to relate loop amplitudes directly to 
string theory in twistor space. Indeed, in the presently known twistor string theories, 
conformal supergravity does not decouple from the gauge theory, which implies that also 
supergravity fields can propagate in loops and the twistor/gauge theory correspondence 
is spoiled at the loop level [19]. 

However, even without the knowledge of the correct twistor string theory description, 
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one can investigate the calculation of loop amplitudes, using recent results. Firstly, the 
recently obtained formulae for tree amplitudes can be fed into the known field theoretical 
unitarity-based methods; and secondly, one can attempt, using twistor intuition, to find 
new field theory prescriptions to compute scattering amplitudes. An example of the latter 
is the off-shell prescription of CSW [6], that uses MHV amplitudes continued off-shell as 
vertices for constructing more general non-MHV amplitudes at the tree level. It is natural 
to try to apply this to the computation of loop diagrams, in particular to the simplest case 
of the supersymmetric one-loop MHV amplitudes, where one can link to known results. 
This will be the subject of this paper. 

Specifically, we will find that twistor-inspired techniques, when applied to one-loop 
diagrams, naturally generate dispersion integrals, in terms of which the amplitudes are 
described. We find that our expression for the MHV scattering amplitudes is in exact 
agreement with the result for this class of amplitudes obtained previously by Bern, Dixon, 
Dunbar and Kosower (BDDK) [29] using the cut-constructibility approach. 

The evaluation of the dispersion integrals proves, rather surprisingly, to be tractable, 
and provides a simpler formulation of the amplitudes in terms of a new form for the 
so-called "two easy masses" box function - see (|7.1)l . We expect this to be an important 
ingredient in further work, both in the calculation of other loop amplitudes, and in the 
development of the twistor space picture. The appearance of dispersion integrals in this 
construction is rather striking. Does twistor string theory presage the return of the 
analytic S-matrix? 

The rest of the paper is organised as follows. In Section 2 we briefly review the 
unitarity-based cut-constructibility approach used by BDDK in [29] in order to derive 
formulae for the one-loop MHV amplitudes. In Section 3 we explain how MHV tree am- 
plitudes, continued off-shell using a prescription equivalent to that of [6], lift to effective 
vertices and can be used to build loop diagrams which are made of collections of MHV ver- 
tices connected by scalar propagators. Loop amplitudes are then obtained by summing 
over appropriate MHV Feynman diagrams. In Section 4 we calculate supersymmetric 
MHV amplitudes at one loop using the procedure outlined in Section 3. Specifically, 
we will combine the off-shell MHV vertices into one-loop diagrams, and show that this 
approach yields dispersion integrals 2 which reproduce the known results for MHV one- 
loop amplitudes (expressed in terms of scalar box integrals). The formal proof of this 
is presented in Section 4. Section 5 is devoted to the explicit calculations of the corre- 
sponding dispersion integrals, for the case of the n-particle MHV scattering amplitude. 
We find a new formulation of the box functions appearing in the amplitude, and show 
that this agrees with the known expressions for the MHV amplitudes at one loop, using 
a nine dilogarithm identity proved in an Appendix. In Section 6 we briefly comment on 
vanishing one-loop amplitudes. Finally, we present our conclusions in Section 7. 

2 For a review of dispersion relations, see, for example, [39]. 



2 



2 One-loop MHV amplitudes from unitarity 



In this section we will briefly review the derivation of BDDK for the one-loop MHV 
M = 4 super Yang- Mills amplitudes from unitarity constraints, given in [29]. We refer 
the reader to this reference for more details. We suppress constant factors connected with 
dimensional regularisation where they are not essential to the discussion. 

The full one-loop n-point MHV amplitudes in M = 4 super Yang-Mills are proportional 
to the tree level amplitudes, 

<r 4MHV = 4TK?- (2.i) 

The function V% is given in terms of scalar box functions F by 

m— 1 n n 

* ' lm 

j 



^2m+l — Fn:r;i + ^n:i 

r=2 i=l i=l 

m—2 n n n/2 
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r=2 i=l i=l i=l 



or, more compactly 3 [25], 



[f]-i 



V " = EE( 1 -2H ! S?- (2 - 3) 

i=l r=l 

The basic scalar box integral J4 is defined by 

f d 4 ~ 2t p 1 

J 4 = -<47T) 2 ^ y ^-)4=* p2 (p _^ l) 2 (p _^ 1 _^ 2) 2 (p + ^ 4) 2 > ( 2 - 4 ) 

where dimensional regularisation is used in order to take care of infrared divergences. 
The relevant integrals arising in the one-loop MHV diagrams are related to J4 for certain 
choices of momenta K, and are denoted by I\™ and If.™f. These are given in terms of 
the F functions in ()2.2)1 by 

o rplrn 

jlm n:i 

L i-3 L i-2 

o p2m e 

T2me nvni ,„ ~\ 

^ Ar+lUr+1] _ ,[r],[n-r-2\ ' \^°) 
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with 

tf = (k + ■ ■ - + fci+r-i) 2 (2.6) 



3 In the next equation it is understood that F%™? = F^™~ (see Eqs. (I.5a)-(I.5b) of [29]). 
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(the ki are the external momenta). The F functions depend on the variables tf for 
certain values of i,r, and involve logarithms and dilogarithms. Appendix I of [29] gives 
the explicit expressions. 

It is also convenient to introduce the general scalar box function [25] 
1 



F(p,q,P,Q) 



+ U 2 (l- 
+ U 2 (l- 

- u 2 (i- 
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P 2 \ f P 2 
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Q 2 
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+ Li 2 (l - 
+ Li 2 (l - 



(P + q) 2 

Q 2 
(Q + p) 2 



log 2 



(P + pf 



(P + p) 2 (P + q) 2 J 2 & \(P + q) 2 
where P + p + Q + q = 0. Introducing the convenient variables 



(P + pf 



t :-- 



(P + qf 



[2.7) 



(2.8) 



we can rewrite (|2.7J) as 
F(s,t,P 2 ,Q 2 ) := 



p2 p2 

LU)- T )+u.(i- T 

P 2 2 \ 1 



Li 2 (l-^ + L H l- t 



St 



(2-9) 



The relation to the functions F 2 ™£ is obtained by setting p = Pi-i, q = Pi+ r , and P = 
Pi + ■ ■ ■ + Pi+r-i- In the following sections we will make use of the discontinuities of the 
function F across its branch cuts. A derivation of these discontinuities up to O(e ) can be 
found in Appendix A and, up to and including 0(e) terms, in Appendix B, where these 
discontinuities are evaluated from phase space integrals. 

The one-loop MHV amplitudes were constructed in [29] from tree diagrams using 
cuts. A given cut results in singularities in the relevant momentum channels, and from 
considering all possible cuts one can construct the full set of possible singularities. From 
this and unitarity one can deduce the amplitude as given in (j2.1|) . More explicitly, consider 
a cut one-loop MHV diagram where the cut separates the external momenta k mi and 
k mi -i, and k m2 and k m2+ i (i.e. the set of external momenta k mi , k mi+ i, k m2 lie to 
the left of the cut, and the set fc m2 +i, k m2 +2, •••> k mx -\ lie to the right, with momenta 
labelled clockwise and outgoing). This separates the diagram into two MHV tree diagrams 
connected only by two momenta l\ and l 2 flowing across the cut, with 



h = h + Pl 



(2.10) 
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where Pl = Y^T 2 ^ ki is the sum of the external momenta on the left of the cut. The 
momenta l%,l 2 are taken to be null. It is important to note that the resulting integrals 
are not equal to the corresponding Feynman integrals where l\ and l 2 would be left off 
shell; however, the discontinuities in the channel under consideration are identical and 
this gives enough information to determine the full amplitude uniquely. 

It is immediate to see that the integrand which arises from the cut diagram involves 
the function [29] 

^ (mi - lm^jkh) (m 2 m 2 + l){hl 2 ) ^ ^ 

(mi - lh){—hmi) (m 2 l 2 ){-l 2 m 2 + 1) 

With the help of the Schouten identity we are able to rewrite 

(m x -1 m x ) (Mi) = (mi-lZi)(Z 2 mi) + (mi-lZ 2 )(miZi) , (2.12) 
(m 2 m 2 + 1) (£1 l 2 ) = (m 2 £ 2 )(£im 2 + 1) + {m 2 h){m 2 + ll 2 ) , 

and hence 

71 = lZ(mi,m 2 + 1) + 1Z{mi — 1, m 2 ) — 7£(mi, m 2 ) — 1Z{mi — 1, m 2 + 1) , (2.13) 
where we define TZ(i,j) to be the homogeneous function of the spinors l\ and l 2 given by 

1 Jj • (</i> o'/ 2 ) ■ 1 j 

We now rewrite the integrand in terms of the scalar functions appearing in the bubble, 
triangle and box integrals. Firstly, we notice that 4 



mi) - { MM£4M - T \v ,.„„ , , 2 j (2.i5) 



i(l+ 7 5 )/iz/ 2 j 

{ik)[ki]{3k)[kj] (h-i) 2 (k+j) 
2[(l l i)(l 2 j) + (l 1 j)(l 2 i)-(l l l 2 )(ij)} 



Of course, one could equally write {l\ i) [h i] = 2{li) = (li + i) 2 = —(li — i) 2 for example. 
The correct choice of the signs in the denominator in 1)2.15)1 is made according to the 
momentum flow. In ()2.15)) and in the following paragraph we have written the signs 
which are appropriate for i = mi, j = m 2 . The other possible cases, corresponding to the 
possible arguments of the TZ functions on the right hand side of (j2.13j) . are treated in a 
similar manner. 



The next step consists in using ()2.10|) to rewrite the second and last terms in the nu- 
merator of the last expression in ()2.15|) . The cancellation of bubble and triangle integrals 

4 In the following formula we will omit a term proportional to an e-tensor contracted with four momenta, 
since it vanishes upon integration. 
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then takes place upon summing over the four terms in (|2.13j) . One may anticipate this 
cancellation by defining an effective function TZ eS (i,j) - this is defined to be 7Z{ij) minus 
the terms which cancel upon the summation involved in (j2.13|) . One finds that 

WW) = ~ 2(I if P ^lf 3) ■ <«<>> 
and one may use this to define a function TZ eG in an analogous way to (j2.13|) . 



Now note the identity 

4{Pi)(Pj) - 2P\ij) = (P + i)\P + jf - P\P + i + jf , (2.17) 

for any momentum P, where i,j refer to null momenta. In other words, the left hand side 
of ()2 . 1 7|) is invariant under 

P^P + ai + bj, (2.18) 

where a and b are arbitrary numbers. Using this, one can show that the momentum 
dependence from the numerator in ()2.16|) precisely cancels that from the denominator in 
the definition of the functions / in (|2.5jl in each of the four cases which arise in TZ cS . 
One is finally left with the result that the cut diagram considered gives rise to the sum 
of four discontinuities (in the same channel) of four different F functions (three when a 
cut leaves only two legs on one side), one for each of the cut box diagrams which arise . 
The knowledge of all cuts, together with the fact that due to general arguments [29] this 
class of amplitudes is given by a sum of scalar box integrals F, is sufficient to fix the full 
amplitude uniquely, and it is simply given by ()2.2|) as a sum over all scalar box functions 
F with all coefficients equal to 1. 



3 Off-shell derivation of the one-loop amplitudes 



The procedure summarized above for calculating the one-loop MHV amplitude involved 
the study of all possible cuts. The calculation of cuts involves an integration over the 
Lorentz invariant phase space (LIPS) measure, which puts the momenta l\, I2 crossing the 
cut on-shell. This fact allows one to insert the tree MHV amplitudes directly into the 
integrals of the cuts, and was also crucial in various algebraic manipulations which led to 
cancellations and the simplicity of the final expression. 

In [6] a prescription for taking momenta of external lines in MHV-amplitudes off shell 
was given, and used to combine tree MHV diagrams in such a way as to produce non- 
MHV tree amplitudes. In this section we will define an off-shell prescription that can be 
applied to one- loop diagrams (and we expect also to higher- loop diagrams). In Sections 
4 and 5 we will show that combining MHV vertices into one-loop diagrams using this 
off-shell prescription precisely yields the MHV results given earlier. 
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Consider an off-shell (loop) vector L. It can be decomposed as [21,23] 



L — I + zrj , (3.1) 

where I 2 = 0, and rj is a fixed (and arbitrary) null vector, rj 2 = 0; z is a real number. 
Equation ()3.1|) determines Z clS cl function of L to be 

(3.2) 



" 2(L V ) • 

Notice that L 2 = z(l + 77) 2 ; since the null vector / has a non- negative energy component 
and 77 is null as well, it follows that (/ + rf) 2 has a definite sign unless the two vectors are 
proportional. This implies that the sign of L 2 is directly related to the sign of z. 5 

We can write the null vectors / and 77 in terms of spinors as U = l a la, Vaa = VaVai 
from which it follows that 6 

[Ifj] 
n a T ■ 

La (l V ) ' 1 j 

These equations coincide with the CSW prescription [6] for determining the spinor vari- 
ables / and / associated with the off-shell (i.e. non-null) four- vector L defined in (j3.1j) . 
The denominators on the right hand sides of ()3.3j) and (J3.4)) will be irrelevant for our 
applications, since the expressions we will be dealing with are homogeneous in the spinor 
variables 77; we will discard them. 

In order to calculate loop diagrams, we need to re-express the integration measure 
d A L (which appears in the expression of loop integrals) in terms of the new variables I 
and z introduced previously. For our purposes it is useful to consider the product of the 
momentum measure and a massless scalar propagator. After a short calculation, one finds 
that 

d L . . - , ,^ dz , . 

_ = d M(l) - , (3.5) 

where we have introduced the Nair measure [40] 

dM(l) := (/ dl)d 2 l - [I dl] d 2 l . (3.6) 

Importantly, the product of the measure factor with a scalar propagator d A L/L 2 of ()3.5|) 
is independent of the reference vector 77. The expression f!3.5|) will be central for our 
construction of loop diagrams. 

5 If I is a null vector, l a a = l a la, then its energy component is Iq = (l/2)(Zi^ + > 0, as in 
Minkowski space we identify la = (la)* ■ 

6 We define the spinor inner products as (A/x) := e^A"/^, [A p\ :— e^A™/)^ 3 . 
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Note that the Lorentz invariant phase space measure for a massless particle can be 
expressed precisely in terms of the Nair measure: 



where, as before, we write the null vector I as l a a = lja> an d in Minkowski space we 
identify 1 = 1*. 

We conclude this section by summarising the strategy that we will follow in order to 
evaluate a generic loop diagram. 

The first step consists in building M HV Feynman diagrams out of MHV vertices. Each 
spinor variable of an MHV vertex corresponding to an internal line is taken off-shell using 
the prescription of CSW outlined in the beginning of this section. Internal lines are then 
connected by scalar off-shell propagators which connect particles of the same spin but 
opposite helicity. Note that each MHV vertex should be multiplied by an appropriate 
delta function for momentum conservation. In the next step we express all the loop 
integration momenta as in (|3.1|) and use the integration measure ()3.5|) which already 
incorporates the scalar off-shell propagators. Finally, one has to sum over all independent 
diagrams obtained in this fashion for a fixed ordering of external helicity states. 

In Section 4 we carry out this programme for the particular case of MHV amplitudes 
at one loop. One of the features of this procedure is the fact that the integration measure 
is naturally expressed as the product of two terms: 

1. a Lorentz- invariant phase space measure, and 

2. an integration over the z- variables introduced according to ()3.1|) (one for each loop 
momentum) . 

The phase-space measure (a two-body phase space measure, at one loop) will be appro- 
priately continued to D = 4 — 2e dimensions in order to deal with potential infrared 
divergences. 



4 Supersymmetric MHV amplitudes at one loop from 
MHV vertices 



We will now write down the expression for the one-loop supersymmetric MHV amplitude 
with n external legs. In performing our analysis, we will make use of the supersymmetric 
formulation introduced in [40], which generalises the usual MHV amplitudes. In this setup, 




(3.7) 
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to each particle one associates the usual commuting spinors X a , Xa (in terms of which the 
momentum of the i-th particle is p aa = A^A^,), as well as anticommuting variables r} A , 
where A is an index of the anti-fundamental representation of SU(4). The supersymmetric 
amplitude can then be expanded in powers of the various Af = 4 superspace coordinates 
rf A , and each term of this expansion corresponds to a particular scattering amplitude in 
M = 4 SYM. A term containing p powers of rf A corresponds to a scattering process where 
the i-th particle has helicity hi = 1 — p/2. 

With this in mind, the expression for a supersymmetric n-valent MHV vertex is 

V(l, . . • , n) = i(2^ ( £ VA*) (£ A y) f[ — ±— , (4.1) 

i=i i=i i=i ' ' 

where the spinors associated to off-shell legs are chosen according to the prescription 
(or, equivalently, ({SHU) an< ^ 




Figure 1: One-loop MHV Feynman diagram, computed in (|4.2|) using MHV amplitudes 
as interaction vertices, with the CSW off-shell prescription. The scattering amplitudes 
with the desired helicities for the external particles are then obtained by expanding the 
supersymmetric scattering amplitude in powers of the J\f = 4 coordinates rf. 

We can now write down the amplitude corresponding to the prototypical Feynman 
diagram we will consider, depicted in Figure 1. This is a diagram where the set of ordered 
external momenta ki, . . . , k n is separated into two sets k mi , . . . , k m2 (on the left), and 
k m2+ i, . . . , k mi -i (on the right). Of course, the full one-loop MHV amplitude is obtained 
after summing over all the possible choices of these ordered sets. We will come back to 
this point later, as it is instructive to first discuss the structure of the generic one-loop 
diagram we want to compute. 

The expression for the Feynman diagram in Figure 1 reads 
A = i{2^\P L + P R ) ( d \ L 2 ld \ L 2 2 S^(L 2 -L 1 + P L ) [dW l dW 2 A L A R , (4.2) 
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where Pl = YlT= mi anc ^ Pr = J^^ml+i are the sum °f the momenta on the left 
and on the right of the diagram, respectively, and the tree amplitudes Al and Ar are 
given by 



m 2 -l 



1 " 1 

(m 2 (2)(-(imi)((2 -ti> .f- 1 - (i« + l> 



i=mi 

mi —2 



= ^'< 9 *> ,„._,,.w-,.„.xiw,.-u n 



(mi - Hi)(-( 2 m 2 + l)(d - i 2 > (jj + 1) 



Here 



(& L )Aa = + V%a ~ V l lha , (4-4) 

ieL 

Notice that (|4.2j) is formally written in four dimensions, but should more correctly be 
analytically continued to D = 4 — 2e dimensions, due to the presence of infrared diver- 
gences. The effect of this is that the Lorentz invariant phase space measure in 1)4.10)1 
below should be taken to be in D = 4 — 2e dimensions. For brevity we will still keep the 
four-dimensional notation. 

Let us stress some important facts about ()4.2j) . 



1. We have introduced off-shell vectors L\ and L 2 as in ljl . i.e. 

Li-a,a licJ-ia ~\~ %i "^]cJ]a j i 1,2. (4.5) 

These only appear in the measure factor d 4: Lxd 4: L2/ \L\L\) , and in the delta function. 
Using (|4.5J) . we will rewrite the argument of the delta function as 

L 2 - L x + P L = l 2 - h + P L , Z , (4.6) 

where we have defined 

P L;Z ■= Pl-z V , (4.7) 

where 

z := z x - z 2 ■ (4.8) 
Note that we use the same 77 for both the momenta Li, % = 1,2. 

2. In the remaining part of the integrand we use, as spinors associated to the off-shell 
(loop) legs with momenta Lx and L 2 , precisely the spinors l\ a and l 2a of (|4.5j) ; this 
is the essence of the CSW prescription. 
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3. We observe that 

^^5^(L 2 -L 1 + P L ) = ^ d -^ dN{h)dN(h) t {A Xh-h+PL, 

L 1 L 2 Z\ Z 2 

dz\ dz 2 



-4 

Z\ Z2 



d%8^\l{) d%5^(l 2 2 ) ^\l 2 -h + P L ., z )\ , (4.9) 

where we used (j3.5j) and ()3.7|) . Note that the expression 

dUPS(l 2 ,-h;P L . z ) := d%6^\l 2 1 ) d%5 (+ \l 2 2 ) 5 (4 \l 2 -h + P L ; Z ) , (4.10) 

which appears in (J4.9j) . is nothing but the two-particle Lorentz invariant phase space 
measure. This will be crucial in the following. We can therefore rewrite ()4.9|) as 

^^8^(L 2 -L 1 + P L ) = -4^^dLlPS(l 2 ,-k;P L;z ). (4.11) 

L l L 2 Z\ Z 2 

To deal with infrared divergences, we will dimensionally regularise the LIPS measure 
appearing in 1)4.11)1 to D = 4 — 2e dimensions. 

Now we return to the evaluation of (J4.2j) . We start off by integrating out the fermionic 
loop variables. This is easily accomplished by first writing 

5 (s \Q L )S (8 \<d R ) = Jd 8 6d s 6' e w «^n e < A {Qn)\ ; (4 . 12 ) 
where and 0# are given in ()4.4|) . Then 

"dW'dW 2 e^f-^)(*'~^ = (1,1^5^(9' - 6) . (4.13) 
In this way, ()4.2|) is recast as 

A = Aree " C , (4.14) 

where the tree-level amplitude Aree is given by 

n n ^ 

Aree = ^(27r)W\P L + P R ) 8® ( £ VaK) II TTTTlT ' ^ 

and the integral £ is 

[d A L 1 d i L 2 {i) . (mi - lffli)(Mi) (m 2 m 2 + l)<Zi / 2 ) , , 

C ~ J^T^l 6 (L2 " Li + Pl) K-l^Km!) (m 2 / 2 )H 2 m 2 + l) " (4 ' 16) 

We recognize in the integrand of 1)4.16)1 the function 7?. defined in (12.11)) . 

We will now proceed to the evaluation of ()4.16)) . Following steps similar to those of 
(J2.11)) - ()2.16j) . we decompose 7Z as in 1)2.13)1 and we are left with a sum of four terms. Let 
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us focus on the term —TZ(mi, m 2 ) in (|2.13|) ; the other cases are treated analogously. This 
term generates a contribution to the amplitude which is proportional to 

fd'L.d'L,^ N(P Z ) 

that is 

1 = / V T ^ ^ " + PiJ ft^WT^ ■ (418) 

where the numerator N(P Z ) is defined by 

iV(P,) := -2(P z -m 1 )(P z -m 2 ) + P*(m 1 -m 2 ) (4.19) 

= -2(P L;2 ■ m x ){P L . z ■ m 2 ) + Pl^mx ■ m 2 ) , 

and 

m 2 — 1 

P 2 := ^ Pi - ^ 77 , (4-20) 

i=m\+l 

Pl,z = Pl ~ zrj = ^ Pi ~ zr l = P z + rrix +m 2 . (4.21) 

i=mi 

Pi = XiXi is the momentum of the z-th particle. 

In the calculation reviewed in Section 2, where the MHV one-loop amplitude is recon- 
structed from the cuts, one was led to integrands involving the function defined in ()2.16|) . 
More precisely, the numerator appearing in the corresponding box integral is the same as 
that in (J4.19)) . but evaluated at z = 0. Equivalently, our numerator 1)4.19)1 is related to 
that in the cut-constructibility picture simply by the shift 

Pl — > Pl,z = Pl-z V . (4.22) 

Next, we change variables from (zi, z 2 ) to (z, z'), where z' = z± + z 2 , so that 

1 - 2 1 { * + %%-,) M[h) si ' ){h ~ l ' + p -» (L-^y+m,)- ' (4 ' 23) 

and perform the integration over z' . Equation ()4.23)) turns into 



/a 7 
— dH{k) dN(l 2 ) 6 {4 \l 2 -h + P L , Z ) — 
z (li 



N{P L ;z 



mx) 2 (l 2 + m 2 y 



-m-r> I— dUPS(l 2 , -h; P L , Z ) N [ 2 "; z) M V2 ■ (4.24) 

' z {Ix - mi) 2 (l 2 + m 2 y 
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Figure 2: The possible cuts of a box diagram, corresponding to the function F defined 
2.9)1 . In the notation of that equation, i and j correspond to the momenta p and q, 



in 



respectively; and P := pi + \ + • • • + Pj-i, Q '■= Pj+i + • • • + Pi-i correspond to the two 
groups of momenta on the upper left and lower right corner, respectively. The vertical 
(horizontal) cuts correspond to the s-channel (t-channel) cuts respectively, and the upper 
left (lower right) corner cuts to the P 2 -channel (Q 2 -channel) cuts respectively. 



In the last step we have used the two-particle Lorentz-invariant phase space measure from 

gU3D. 

Let us briefly comment on the appearance of the LIPS measure in (|4.24|) and its con- 
sequences. This is the point where the off-shell calculation presented in this paper makes 
contact with the approach of BDDK, where one-loop amplitudes are reconstructed from 
the evaluation of the cuts in the various channels. Indeed, the LIPS measure is precisely 
what is required by Cutkosky's cutting rules [41] to compute the discontinuity of a Feyn- 
man diagram across the branch cuts. Which discontinuity is evaluated is determined by 
the argument in the delta function appearing in the LIPS measure; in the case of (|4.24j) . 
this is P L - Z . 

Specifically, the phase space integral appearing in ()4.24|) is computing a particular 
discontinuity of a box diagram. The generic box diagram, corresponding to the function 
F defined in ()2.9|) is depicted in Figure 2. The phase space integral appearing in (|4.24p . 
which in turn was generated by the term lZ(mx,m2), corresponds to the particular (cut) 
box diagram with i = mi and j = mj, and where the momentum flowing in the cut is 
equal to J2ml V% — zr l = ~ (Q + zr i)- ft corresponds to the lower right corner cut in Figure 
2, where however the momentum Q has been shifted to Q z := Q + zrj (and consequently 
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P — > P z := P — zr], so that momentum conservation reads P z + Q z +P+Q = -P + Q+p+g = 
0). Hence we conclude that the phase space integral appearing in the last line of (|4.24|) 
is computing the (^-discontinuity of the box function F(s z , t z , P 2 , Q 2 ), where 

P z := P-zri, Q z := Q + zrj , (4.25) 

and 

s z := {P z +pf , t z := (P, + g) 2 . (4.26) 

It is now important to remember that the full result for the one-loop calculation of the 
MHV amplitudes includes a sum over various MHV Feynman diagrams. To generate all 
the possible diagrams, we rotate the external momenta in Figure 1, so that the position 
mi in Figure 1 will be taken by all possible external momenta Pi, ■ ■ ■ ,p n - Moreover, we 
should also vary the number of momenta on the left (and hence on the right) of the loop 
by varying the momentum appearing in the 777,2 position. For each diagram, we will follow 
the same steps as discussed above; in particular we will apply the Schouten identity twice, 
so as to generate four terms from each diagram, as in (j2.13|) . 

Thus we note the following two key consequences: 

1. Firstly, in this way we will produce, for each fixed box function F, exactly four 
phase space integrals, one for each of all the possible cuts of the function; and 

2. secondly, the functions F generated by summing over all MHV Feynman diagrams 
are precisely those of the double sum of ()2.3j) . where the sum over i corresponds 
to different choices for the momentum in the rri\ position, and the sum over r 
corresponds to varying the number of momenta between the legs mi and 777 2 in 
Figure 1. Each function F appears precisely once (in all possible cuts). 

Hence, in order to prove that our procedure correctly generates the known MHV one-loop 
amplitudes ()2.2|) . it suffices to focus on a single function F. 

Keeping in mind the previous considerations, we now come back to ()4.24|) . where the 
Lorentz invariant phase space integral represents one of the four possible cuts of the F 
function we wish to reproduce. These Lorentz invariant phase space integrals correspond 
to the discontinuities of the box function F(s z ,t z , P 2 ,Q 2 Z ) in the four "channels" s z , t z , 
P 2 , Q 2 . Each of the four cut-box functions is then integrated over z. 

Now we discuss how the final z-integral of the sum of the four cuts of F(s z , t z , P 2 , Q 2 ) 
reproduces the function F(s,t, P 2 ,Q 2 ), hence proving that the procedure discussed here 
correctly generates the MHV amplitudes at one loop. 

Consider for example the s z -channel, for which the variable Pl ;z appearing in ()4.24j) 
is equal to P + p — zrj. We define 

s ' ■= Sz = (p + p - zr] ) 2 = s - 2zr)- {P + p) , (4.27) 
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where the kinematical invariant s is given by s = (P + p) 2 (see (I2.8jl ). It then follows that 
ds' = —2dz(r]Pi J ), and therefore 



dz ds' . , s 

- = • 4.28 

z s' — s 

This last step allows us to recast each of the four terms as a dispersion integral; the 
dispersion integrals reconstruct the function F from its discontinuities. The sum over 
Feynman diagrams built out of MHV vertices corresponds to summing dispersion integrals 
for all possible channels for each function F appearing in the sum ()2.2|) ; and then summing 
over all the functions appearing in (j2.2j) . In this way one reconstructs the full MHV 
amplitude (|2.1j) at one loop. 



We will now show in more detail how the function F(s, t, P 2 , Q 2 ) is reconstructed by 
summing over dispersion integrals. We denote by T the result of the ^-integration of the 
sum of the four phase space integrals (each of which is the discontinuity of F in one of its 
variables), that is 

T := Ids'^L + + fdP 2 '^- + fdQ 2 '^- . (4.29) 

J S'-S J t'-t J P 2'_p2 J ^ Q2'_Q2 ^ > 

As is customary, for a function f(z) which is analytic in the z-plane except for a cut on 
the real axis at z > zq G R, we define the discontinuity of f(z) across the cut as 

A z f := f( z + ie)-f( z -ie) , (4.30) 

where e — >• + and z > z . 



Now we want to show that T is actually proportional to the function F itself. More 
precisely, we will show that the discontinuities of the function T in s, t, P 2 and Q 2 
are the same as those of the function F, modulo a proportionality coefficient. One can 
further argue that, in supersymmetric theories, scattering amplitudes can be entirely 
reconstructed from the knowledge of the discontinuities [30] - no ambiguities related to 
the presence of rational functions occur in this process - to conclude that T must actually 
be proportional to F. 

To demonstrate the equality of the discontinuities, we notice that discontinuities on 
the right hand side of ()4.29j) come potentially from two sources: 



1. Discontinuities due to the integral. Those are calculated by simply remembering 
that 

1 



x ± is 

where P stands for the principal value prescription. 



P — =F m5(x) , (4.31) 
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2. Discontinuities in the functions A S >F, A t /F, A p2 'F, A Q2 >F. 

However, we can immediately rule out discontinuities of type 2. The explicit expressions 
of the discontinuities A S >F, A t >F, A p2 >F, A Q2 >F are listed in the Appendix, and are 
expressed in terms of logarithms. The key point is that, for each integral, the corre- 
sponding arguments of the logarithms appearing are always positive; no discontinuities 
can therefore be generated. We are thus left with the discontinuities of the type 1. Using 

— — = -2m5{x) , (4.32) 

x + is x — is 

we immediately get 

2vriA s F, (4.33) 
2-ni A t F , 
2ixi A P 2 F , 
2ixi A Q 2 F . 

We have found that the functions T and F have precisely the same cuts in all channels 

(s, t, P 2 , Q 2 ). 7 By appealing to the cut-constructibility assumption, we conclude that 

T = —2ni F . (4.34) 

Finally, notice that by picking the 5-function contribution in (|4.31|) . we lose any depen- 
dence of the discontinuities on the arbitrary vector rj which was used in (13.1)1 to write the 
expressions for the non-null loop momenta L\ and L>2- 

In the above, we have presented a general argument showing that if one combines 
MHV vertices with a suitable off-shell prescription, one generates the expression f)4.29j) . 
giving the box function T in terms a sum of dispersion integrals. The full n-particle 
amplitude is then given as a sum over box functions, as specified in Section 4. 

One may wonder if the presence of infrared singularities in the integrals affects this 
general argument. In the next section we will present the calculation for the n-gluon 
scattering amplitudes at one loop, showing how the integrals of the form 1)4.24)1 give rise 
to those in (|4.29j) . how the singularities in these integrals are dealt with, and finally 
evaluating explicitly the sum of dispersion integrals appearing in (j4.29j) . 



A S T = - 

A t T = - 

A p2 J" = - 

AniF = ~ 



5 The n-particle scattering amplitude at one loop 

In this section we will consider the general MHV, one-loop, n-gluon amplitude, and show 
that the calculation of the dispersion integrals ()4.29j) directly yields the box function (J2.7j) . 

7 We observe that s + t + u = P 2 + Q 2 , where u — (p + q) 2 . 
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Summing over the MHV Feynman diagrams (Figure 1) obtained by gluing MHV vertices 
at one loop using the procedure detailed in the sections above then gives the sum over 
box functions ()2.3|) which represents the full MHV, one-loop, n-gluon amplitude. 



As explained in the previous section, we can focus on a single box diagram. The goal of 
this section is to show that the formal expression (|4.29j) really reproduces F(s, t, P 2 , Q 2 ). 



In Appendix B we have computed the phase space integrals which give the discon- 
tinuities of the box diagram represented in Figure 2 in the channels s, t, P 2 , Q 2 . The 
main result is displayed in (|B.17j) . It then follows that the dispersion integral in, say, the 
s-channel is given by 



I(s) 



ds' 



s — s' 



as'Y 



1 + 

\1 — as'/ 



where 



P 2 + Q 2 - s - t 
P 2 Q 2 - st 



a 



P 2 Q 2 



st 



(5.1) 



(5.2) 



A few comments are in order: 



1. In writing (|5.1j) we have used the form (|B.17|) of the discontinuities of the box 
diagram, where we must also keep the 0(e) dilogarithm term in brackets in (|B.17|) . 
This is because the dispersion integrals potentially give rise to 1/e divergences, and 
hence the discontinuities have to evaluated up to order 0(e) if the amplitude is to 
be calculated up to order 0(e°). 

2. In the integrand of ()5.1|) we have omitted a numerical factor that depends on the 
dimensional regularisation parameter e. This factor, which is explicitly written in 
(|B.17J) . is irrelevant for our discussion. 

3. We have chosen the reference vector rj to be equal to p. We know from the previous 
section that our final result will be independent of r], at least for the cut-constructible 
part. However, what we suggest here requires a stronger gauge invariance, namely 
that rj can be chosen separately for every box function, i.e. rj is fixed for the four 
integrals of the type ()5.1|) contributing to a particular box, but can be chosen inde- 
pendently for any box. 8 Indeed it turns out that the rj dependence in the integral 
()5.1|) cancels when combined with the integrals for the other three channels. We 
have proven this fact for the terms, which are singular and finite in e, by numerical 
integration for various arbitrary Euclidean points, i.e. s, t, P 2 and Q 2 all negative, 

8 A related issue was discussed in [6] to establish independence of tree-level amplitudes from the 
reference spinor rf needed in the definition of the MHV vertices. There gauge independence is only 
recovered after summing over all MHV diagrams, where rj a has to be kept fixed for all diagrams, and not 
only for a subset. 
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and varying rj randomly. 9 The particular choice rj = p (or r\ — q) has some ad- 
vantages. One of them is that for this choice, it turns out that the quantity N 
introduced in (jB.21|) becomes extremely simple and independent of z. Indeed, one 
has 



s z t z - P z z Qi = st- P l Q l + Az 
from which it follows that 



(pq) (rjP) - (Pq) (rjp) - (Pp) (rjq) + 4z 2 (r)p) (rjq) , 

(5.3) 



s z t z - P'Qi = st - P Z Q Z , for rj = p or q . (5.4) 

In the following we will choose 77 = p. For this choice, the quantity a defined in (|B.18|) 
becomes a constant, and its expression simplifies to ()5.2|) . 

The next important observation in our calculation is that the combination of dispersion 
integrals from ()4.29|) is equal to 

T = I(s) + I(t) - I(P 2 ) - I(Q 2 ) . (5.5) 



Considering the combination I(s) — I(P 2 ), expanding the expression (1 — as') e in 
powers of e, and using Landen's identity 10 



Li2 (l3^) + \^g 2 (l- as') = -Lisas') 



(5.6) 



one finds that 
I(s)-I(P 2 ) - 



s-P 



2 poo 



ds' 



{s-s'){P 2 -s 



-(sT e l+elog(l-as')-e 2 Li 2 (as') . (5.7) 



Now consider the three terms in the integrand above in turn. The first may be directly 
integrated to yield 



1 



Ti = -~ [en csc(e7r)] (s)-* - (-P 



2\-e 



(5.J 



The second term in the integrand of ()5.7|) gives 

ds' 



To 



s - P l 



(s - s')(P 2 - s') 



logfl — as') 



(5.9) 



1 r 



-a) e {aP 2 )- £ n {as)' 6 csc(evr) {aP 2 ) e {asy(H{e,aP 2 ) - H{e,as)) 



+ e(as) e log(l - aP 2 ) - e(aP 2 ) 6 log(l - as 



9 We thank Lance Dixon for prompting us to clarify this point and for sharing his results on an 
(independent) numerical proof. 



10 



This form of Landen's identity applies for as' (1, 00). 
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where we define H(e, z) = 2 -Fi(l, e, 1 + e, 2), with 2 -Fi the hypergeometric function. Now 
note that 



H(e,z) = 2 F 1 (l,e,l + e,z) = (l-z)-* 2 F 1 (e,e,l 
= (l-z)- e (l + e 2 Li 2 (^ I )+0(e : 



z- 1 



(5.10) 



Using this and expanding in powers of e, one finds that 

aP 2 



T 2 = [ire csc(7re)] 



Li? 



aP 2 - 1 



Li 2 



(IS 



+ Ilog 2 (l-aP 2 )-ilog 2 (l-a S ) 
as — 1 / 2 2 



\og(aP 2 ) log(l - aP 2 ) + log(as) log(l - as) 



(5.11) 



Then, noting the dilogarithm identities 11 

T . / X \ 1 

Ll 2 



+ -log 2 (l-x) = -Li 2 (x) 
x — 1 / 2 



-Li 2 (x) - log(x) log(l - x) = Li 2 (l — x) - 



71 



one finds that 



T 2 = [7re csc(7re)] Li 2 (l — aP 2 ) — Li 2 (l — as) 



(5.12) 
(5.13) 



The third term in the integrand of (J5.7)) gives T^{s) — T^P 2 ), where 

ds' 



TJs) = -( 



s — s 



7 ( S ')- e Li 2 (a S ') 



However, 
Ts(s) 



— — S 67T CSC 67T 



^as) £ LerchPhi(as, 2, e) - Li 2 (as) = C?(e) , (5.15) 



(5.14) 



the last equality following from the fact that LerchPhi(x, 2, 0) = Li 2 (a;). Thus the result of 
this part of the integral in ()5.7|) vanishes as one takes e = 0, and hence it can be dropped. 

Collecting the above results ()5.8j) . 1)5. lHj) . ()5.15j) . together with similar equations with 
s and P 2 replaced by t and Q 2 respectively, we conclude that the expression for the box 

11 The first of these relations is again Landen's identity (|5.6(l : the second one is due to Euler. 
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function F, as obtained in our approach, is equal to (in the following expression we drop 
the overall, ubiquitous factor 7recsc(7re)) 

1 



F{s,t,P\Q*) = -- (-a)- + (-t)-< - (-P 2 )- - (-Q 2 ^ 

+ Li 2 (l - aP 2 ) + Li 2 (l - aQ 2 ) - Li 2 (l - as) - Li 2 (l - at) , 

(5.16) 

p 2 + q 2 - s - 1 _ u 

G " P 2 Q 2_ st " P 2Q2_ st ■ ^- U > 

This expression f)5.16j) should now be compared with the previously known form for the 
function F given in ()2.9j) . One then concludes that (|5.1fij) and (|2.9jl coincide whenever 
the following equality is satisfied: 

Li 2 (l - aP 2 ) + Li 2 (l -aQ 2 ) -Li 2 (l- as) - Li 2 (l - at) = ^ log 2 (~ 
Li, ( 1 - — )+ Li 2 (\ - — )+ Li 2 (\ - ^ V Li 2 ( 1 - ^ V Li 2 ( 1 P2( ^ 



where 



s / V i / V s I V t / V st 

(5.18) 

This is a remarkable identity involving nine dilogarithms. We discuss and prove it in 
Appendix C. Here we want to mention that a region in the space of kinematical invariants 
where the identity (|5.18jl holds is when s, t, P 2 and Q 2 are all negative, i.e. the Euclidean 
region. In this region the dispersion integrals can be computed safely and, moreover, none 
of the logarithms and dilogarithms in ()2.9|) have their arguments on their respective cuts. 
The representation ()5.16|) of the function F, which therefore coincides with ()2.9|) in this 
region of the space of kinematical invariants, can then be analytically continued to generic 
values of the invariants. Let us also point out that another interesting region where (|5.18j) 
holds occurs when one of the invariants, say s, is taken to be positive and shifted by a 
small positive imaginary part, s — >• s + ie, e — > + , while all the other invariants are 
negative. 

In order to obtain the one-loop amplitude in a physical region, one has to perform an 
analytic continuation. Usually, this continuation is simply achieved by the replacement 
s — » s + ie, e — > + , where s stands for any one of the kinematical invariants. However, 
it was pointed out in [45] that this procedure can fail if a cut is hit by the logarithms or 
dilogarithms. In particular, the fith logarithm in ()2.9|) is problematic [45] and has to be 
amended by correction terms after analytic continuation, see Eq. (A.5)-(A.6) in [45]. On 
the other hand we have checked that our form of the box function (|5.16j) does not suffer 
from this problem and, hence, is valid in all kinematical regions. 

Summarizing, we saw in previous sections that combining MHV vertices into one- 
loop diagrams yields the dispersion integral representation (|4.29jl for the box functions. 
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Performing the dispersion integrals explicitly, we have shown in this section that one 
obtains the result f)5.16|) . Using the equality (j5.18|) . one sees that the expression for the 
box function (j5.16|) precisely reproduces the known box function ()2.9|) . Therefore we find 
agreement with the known expressions [29] for the MHV amplitudes at one loop. 

We also note that the expression of the function F given in (|5.1fi|) contains only four 
dilogarithms and no logarithms compared to ([2.9)1 . which contains five dilogarithms and 
one logarithm. 



6 Vanishing one-loop amplitudes 

In this section we make some brief remarks on the vanishing of the simplest one-loop 
amplitudes, and their construction from MHV vertices. 

To begin with, consider the one loop amplitude ( ) with two external gluons with 

negative helicity. This amplitude violates helicity conservation and should therefore van- 
ish. Indeed, we can easily see this by an explicit computation of the diagram contributing 
to the process built out of the MHV vertices (continued off-shell as proposed by CSW). 
It is instructive to perform the calculation in a gauge theory with generic matter content 
(hence without making reference to the Nair M = 4 supervertex) . Then, summing over 
allowed helicity assignments on internal legs and over allowed particles which can circulate 
in the loop, one finds that the amplitude is proportional to a factor of 

1 - m + n s , (6.1) 

where rif is the number of Weyl fermions and n s the number of complex scalars in the 
theory. This vanishes for any supersymmetric theory. 

By applying the procedure described in Section 3, it is also easy to show that the 
M = 4 amplitudes (— — • • • — ) vanish at one loop. To this end, recall that for an L-loop 
amplitude with q external gluons with negative helicity, the number V of MHV vertices 
which are necessary to construct it is V = q — 1 + L. At one loop, an amplitude with n 
negative-helicity gluons requires precisely n MHV vertices. A typical diagram contributing 
to the n negative helicity gluons amplitude is depicted in Figure 3 (for n — 5). It is then 
immediate to see that all possible cuts one can make on the loop diagram give rise, on 
both sides of the cut, to tree amplitudes which are zero. Indeed, these tree amplitudes 

must be of the form H • • • — and hence vanish. A similar reasoning works for the 

(+ — • ) amplitude. 

Finally, we remark that assuming that amplitudes should be constructible from cut 
diagrams and MHV tree vertices, the amplitudes (=F +•••+) are then trivially zero at 
one loop since they cannot be made from two or more MHV vertices. 
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Figure 3: One of the MHV diagrams for the ( ) process. The tree-level diagrams 

on both side of the cuts are of the type — • • • — -\- ; and hence vanish. 

7 Conclusions 

In this paper we have shown that supersymmetric MHV amplitudes at one loop can be 
obtained by using the tree level MHV diagrams as vertices, complemented by the CSW 
prescription to take some of the external lines off shell, when these lines form part of 
an internal loop. The Feynman rules for the combination of MHV vertices into one-loop 
diagrams yield a sum over terms which precisely corresponds to the sum over terms in the 
dispersion relation realization of the one-loop amplitude. This appears rather remarkable, 
and is an appealing feature which one might expect to persist in more general calculations. 

In particular, one would naturally seek to extend this work to compute next-to- 
maximal helicity violating amplitudes in M = 4 super Yang-Mills at one loop using MHV 
vertices, which are only known up to n = 6 external legs. Further applications would 
be to compute more general non-MHV amplitudes for an arbitrary number of external 
gluons, as well as the computation of higher-loop gauge theory amplitudes, and the study 
of theories with less supersymmetry. 

Let us also mention that our analysis provides a representation (j5.16|) of the "two easy 
masses" box function F in ()2.9j) expressed in terms of only four dilogarithms and with no 
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logarithms: 



F(s,t,P 2 ,Q 2 ) 



Q 2 )- e 



(7.1) 



+ Li 2 (l - aP 2 ) + Li 2 (l - aQ 2 ) - Li 2 (l 



as) - Li 2 (l 



at) 



where 



a = 



P 2 + Q 2 - s - t 
P 2 Q 2 - st 



P 2 Q 2 - st 



u 



(7.2) 



It seems likely that this new, simpler representation will provide clues as to the correct 
twistor space description of M = 4 super Yang-Mills, and might simplify the form of other 
loop amplitudes. 

Finally, it is intriguing that the one-loop expression (J2.3)) for the MHV amplitudes is 
a sum of box integrals - each with coefficient one. In term of the construction with MHV 
vertices, this comes only after elaborating the initial expression (by multiple applications 
of the Schouten identity), and appropriately collecting cut box integrals in order to recon- 
struct, from the dispersion integrals of the different cuts, the full box function. It would 
be remarkable if one could find a (string?) theory where the box diagrams represented in 
Figure 2 themselves have a more fundamental and direct meaning as "Feynman diagrams" 
emerging from a perturbative description. 
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Appendix A: Discontinuities of the function F: a first 
calculation 



In this Appendix we compute the discontinuities of the function F(s,t, P 2 ,Q 2 ) defined 
in (|2.9J1 considered as a function of the kinematical invariants s, t, and the external 
momenta P 2 and Q 2 . We should immediately say that in this way we will miss (part 
of) the 0(e) terms in the discontinuities of the box diagrams corresponding to the box 
functions (which of course does not include all terms which vanish as e — > 0). As explained 
in Section 5, these 0(e) terms are important and should be included in the evaluation 
of the dispersion integrals. The complete calculation of these terms requires the explicit 
calculation of Lorentz invariant phase-space integrals, which is presented in Appendix B. 
It will nevertheless be instructive to first calculate directly the discontinuities of the F 
function. 

Each discontinuity is computed in the following way: all kinematical invariants are 
taken to be negative except for one, say s, which is positive (i.e. in the physical region); 
the discontinuity in a certain function / of the invariants is then computed as A s / = 
f(s + is) — f(s — ie), with e — > + . 

The discontinuities are then given by: 



A s F 



+ 




(A.l) 



At F 



+ 




(A.2) 



A p2 F 




(A.3) 



1-ni-\P- 



e 



+ 2iri log 



( 



(s-P 2 )(t-P 2 
(st - P 2 Q 2 ) 



) 
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F 



-,\Q 



log(l- 



~ e 2isin(7re 
P 2 Q 2 



2m 



log(l-^) +log( 



i-9l 
t 



St 



(A.4) 



+ 27Tzlog 



>-g 2 )(t-g s 

(st - P 2 Q 2 ) 



We notice that, for each "channel", the corresponding discontinuity is perfectly well de- 
fined - the arguments of the various logarithms are always positive. For example, in the 
s-channel one has s > 0, all the other kinematical invariants (t, P 2 , Q 2 ) being negative; 
and similarly for the other channels. 



Appendix B: Discontinuities of the function F from 
phase-space integrals 



Consider the phase-space integral which occurs in (j4.24j) . 

V= [d D UPS(l 2 ,-h;P L;z ) T N [^: z \ V2 , (B.l) 

where we have introduced dimensional regularisation in dimension D = 4 — 2e [42] in 
order to deal with infrared divergences, with 

rf D LIPS(/ 2 ,-/ i; P L;2 ) := d D l 1 5 (+ \l 2 1 ) d%5M(lj) 6W(h-h + P L . t ) . (B.2) 

The numerator N(P L . Z ) is defined in (j4.19|) (and is a constant with respect to the phase 
space integration). 

The phase-space integral (|B.1J) computes the discontinuity of the box diagram repre- 
sented in Figure 2 in the P\. z channel. In the following we will use the results of [43,44] 
to evaluate (jB.lj) . 

To begin with, we observe that the delta function in (|B.2)1 localizes the integral (jB.lj) 
onto vectors l 2 , h such that l\ — l 2 = Pl ;z - Lorentz transform to the centre of mass frame 
of the vector l\ — l 2 , so that 

h = \pl ;z {-1, -v) , h = \pl, z {-1^) , (B.3) 

and write 12 

v = (sin 6i cos 6 2 , ••• , cos#i) . (B.4) 

12 We refer the reader to Appendix B of [44] for a careful treatment of the dimensional regularisation 
of integrals such as (|B.1|I . 



25 



Using a further spatial rotation we write (^4, B, C are constants and the difference between 
the four vector mi and the component mi will be apparent from the context) 



mi = (mi, 0, 0, mi) , m 2 = (A, B, 0, C) , 
with the mass-shell condition A 2 = B 2 + C 2 . With this parameterisation, 
(h -mi) 2 (l 2 + m 2 ) 2 = -4 (hmi) (l 2 m 2 ) 

2 



(B.5) 



(B.6) 



4 



P, 



L:z 



mi (1 — cos 0i)(A + B sin 9i cos9 2 + C cos#i) . 



A short calculation shows that, after integrating over all angular coordinates except 9i 
and 6*2, the two-body phase space becomes 



cf- 2e LIPS(/ 2 ,-/ i; P L;2 



7T2 



4r(|- e ) 

Using (jE3 and (|B~6l) . we recast (IB~Tll as 



-2c 



d9id9 2 (sm9iy- 2e (sm9 2 )- 2e . (B.7) 



V = A 



7T2 



4r(i-e) 



-2c 



J 



(B.8) 



where A = N(P L . Z )/Pl. z mi, and is the angular integral 
J := I d9i / d9 2 



2t (sineO^rsin^) -26 



1 - cos 0i)(A + C cos 6»i + P sin 6»i cos 6» 2 ) 



;b.9) 



This integral (|B.9|) has been evaluated in [43]; we borrow its result in the form of [44], 
with the result 



J 



2A 



e 2 Li 2 



A-C 
2A 



0(e 



A + C \e J \A + C / 
To express results in Lorentz covariant form, we note that 

N(P L;z ) = -Pl. z (A + C)m u mi-m 2 = mi(A - C) 
from which it follows that 



1 



A + C 

2A 
A + C 

2A 
A-C 



.V(7',, : ) 

(mim 2 )Pl 



1 

A 



N(P L ; Z ) ' 

N(P L ; Z ) 

(m 1 m 2 )P| • 



(B.10) 

(B.H) 

(B.12) 
(B.13) 
(B.14) 
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Putting together the above results we find that 



V 



3 

7T2 



r(i-e) 



Pi 



2A 



A + C 



1 + e 2 Li 2 



2A 



+ 0(t 



(B.15) 



To connect this with (|4.29jl . let m,\ = p, m 2 = q, and note that the quantity (j4.19|) can 
be rewritten as 

(B.16) 



2N(P Z ) = P 2 Q 2 Z - s z t z . 



Using the expressions for 2 A/ (A + C) and (A - C)/2A given in l)B.13|) . ()B.14jl (with 
mi — > P an d ~^ ?)) we can recast (jB.15|) in the form that we have used in calculating 
the n-gluon MHV scattering amplitudes: 



V 



3 

7T2~ 



r(i-«) 



p, 2 



(i - 



where 



(pg) 



1 + e 2 Li 2 



1 " ^l, Z 



+ 0{e 



II 



(B.17) 
(B.18) 



and 



N{P L , Z ) 2N{P L , Z ) ' 
u = (p + q) 2 = P 2 + Q 2 -s-t . (B.19) 

We also notice the following identities, associated to the s, t, P 2 , Q 2 cuts respectively: 

(B.20) 



- (s z 


~P 2 z ){s z 


-Ql) 


= 2(P^p-q- 


n(pH)) , 




-P!)(t z 


-Ql) 


= 2(P 2 ^p.q- 


N(P$ Z )) , 


(s z 


-P 2 z ){t* 


-Pi) 


= 2{P 2 L { pp.q 


-n(pC)) 


{s z 


-Qt)(t z 


-Ql) 







where P^] = P z + p,P z + q,P z ,Q z respectively, for labels a = s,t, P 2 ,Q 2 . Note that 
these identities hold for any choice of the spinor field rj which arises in the definition of 
the off-shell loop momenta. We also have 



N(P J 



('h 

L:z/ 



N(P J 



(t)^ 
L\z) 



N(P t 



N(P J 



L;z > 



\{PlQl 



S Z tz) 



(B.21) 
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Using ()B.20j) and (|B.21J) . we get the following results for the expression of 2 A/ (A + C) in 
(IB.13|) in the s-, t-, P 2 - and Q 2 -channels: 

1 ( P lITp ■ 1 _ {s z -P*){s z -QD 
N(P%) ~ s z t z -P-M ' 



(p£ z ) 2 p-q (t z -p 2 )(t z -Q 



1 



2^ 



N{P® ) s z t z - P*Ql 

(pV) 2 p ■ V {s z -Pj){t z -Pt 



N(P^) s z t z - PIQl 



[pfjfp-q (s z -Ql)(t z -Q 



1 - 



2^ 



N(Pt Q P) s z t z - PIQl 



The expressions on the right-hand sides of ()B.22|) are precisely those appearing in the cuts 
of the box function F, as given in the log terms in equations (A.1)-(A.4) of Appendix A. 
Thus we see that the LIPS integral in ()4.24|) generates the cuts AF in the box functions. 
For each cut variable a, we then change variables from z to a' — {P^ z ) 2 - Then the sum 
of the four integrals in ()4.24|) corresponding to the four cuts yields the expression ()4.29|) 
as discussed earlier. 



Notice that our result for V in (|B.17jl contains also an 0(e) term (the dilogarithm 
term on the right-hand side of (|B.17j) . This term could not be obtained from the direct 
calculation of Appendix A, and, as explained, is important in the calculation of dispersion 
integrals of Section 5. 

For reference, we write down the different values of ((A — C)/2A)~ l (where ((A — 
C)/2A) is the argument of the dilogarithm in (jB.15|) ) in the different channels using 
dEH and (IB~20l) : 





i 


{s z -P 2 z ){s z 


-Qz) 


{P ( lI?P ■ 


q 


s z u 






i 


(t z -p 2 ){t z 


-Ql) 


(p£ z ) 2 p ■ 


q 


t z u 




n(pV: 


i 


{s z -P 2 z ){t z - 


P 2 Z ) 


(C } )v 


q 


p> 




mpffP: 


i 


(s z -Q 2 z )(t z - 


Ql) 


(< } )V 


q 


Qlu 





(B.23) 



The expansion of ((A + C)/2A) e in (jB.15|) in powers of e leads to powers of logarithms 
with argument 2A/{A + C). These arguments are listed in ()B.22j) for the four possible 
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channels. The arguments of the dilogarithm in (|B.15J) for all possible channels are listed 
in ()B.23|) . We notice that, for each channel, whenever the argument of the logarithm is 
positive, the corresponding argument of the dilogarithm is greater than 1. In other words, 
in the kinematic regime where the phase-space integral is evaluated, both functions (log 
and dilog) are continuous functions of their arguments. 13 

Appendix C: Proof of the identity with nine 
dilogarithms 

In the calculation of the n-gluon scattering amplitude presented in Section 5 we have 
encountered the interesting identity (|5.18jl involving nine dilogarithms. The purpose of 
this Appendix is to give a proof of this identity. 

It is instructive to present the proof for three separate CclSGS, clS they rely on different 
dilogarithm identities: 14 

(a) P 2 = Q 2 = 0; this is the case relevant for the four-gluon scattering amplitude; 

(b) P 2 7^ 0, Q 2 = 0; this is the case relevant for the five-gluon scattering amplitude; 
and finally, 

(c) P 2 7^ 0, Q 2 7^ 0; this is the case relevant for the n-gluon scattering amplitudes, with 
n > 6. 



We start off by addressing case (a). In this case, 

a \p2=Q2=o — ~ + 7 > (CI) 
where a is defined in (15.17)1 . Using Li 2 (l) = vr 2 /6, we see that ()5.18j) becomes 

-i*(-£)-i*(-l) = I V (f) + f (c) 

This equation is satisfied thanks to the identity (z ^ (0, 1)): 

Li 2 (z) + Li 2 Q + ^log 2 (-^) + y = . (C.3) 



13 In our conventions, the cut of the logarithm is on the real negative axis; consequently, the dilogarithm 
has a cut on the interval (1, oo). 

14 A list of dilogarithm identities can be found in the relevant section of [46]. 
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Next we discuss case (b). In this case we make use of the relations: 

u 

l-as\ Q2=0 = 1 + - , (C.4) 



l-at| Q2=0 = 1 + 



]-aP 2 L , „ = (l--)(l- P2 



,2 -o v s J V t 



Now the identity (|5.18j) . which we wish to prove, turns into: 



Li 



Li 2 ( 1 — — j-Li„(l- — ) = (C.5) 



Li2 ( 1 + ^) +Li2 ( 1 + ^) + _ log2 ( 
Eq. ()C.5|) is the same as Hill's five dilogarithm identity: 



Li 2 ( Z uO-L l2 ( Z )-Li !(tt O-Li 2 (^^)-Li 2 (^<^)-ilog 2 (i^) = 0, (C.6) 



if we make the following assignments: 



p2 p2 

— , w = 1 - — , (C.7) 
s t 



which imply 



w(z—l) P 2 — t u z(w — l) P 2 — s u 

! + -> "I " = 7 = 1 + T> ( C - 



I — w s s 1 — z t t 

and 

1 — w s 



1-z t 



(C.9) 



Finally, we consider the general case (c). The relevant identity will be Mantel's identity 
involving nine dilogarithms: 



Li O - U ^)M1) +U '{y) +U ^) < ai °> 
+ Li 2 (x) + U 2 (y) - Li 2 (v) - U 2 (w) + l - log 2 ( - -) , 

where 

(l-v)(l-w) = (l-x)(l-y) , (C.ll) 
and all variables v, w, x, y, lie on the real axis between and 1. 

In order to use it, we notice that if we choose 

x = as , y = at , v = aP 2 , w = aQ 2 , (C.12) 



v\ T . (W 
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then (jC.llJ) is satisfied; in particular, 

x (s - P 2 )(s - Q 2 )(t - P 2 )(t - q 2 ) , 
(! - «)(1 - «>) = (1 - x)(l - y) = { - {st _ p2 Q 2 Y ■ ( C - 13 ) 

Now we use Mantel's identity (jC.10|) with the assignments (|C.12[) . to get: 

+ Li 2 (as) + Li 2 (at) - Li 2 (aP 2 ) - Li 2 (aQ 2 ) + - log 2 ( - = . 

To relate the identity (|C.14|) to our identity (j5.18J) . we need a relation which connects 
Li 2 (z) to Li 2 (l — z). This is Euler's identity, 



TV 2 



U 2 (z) = -Li 2 (l - z) - \og(z) log(l - z) + — . (C.15) 

o 

By repeatedly using (jC.15|) in ()C.14|) . we get: 

Li2 1 t _ ^ V Li2 U _ H\ _ Li2 A - ^ - Li 2 f 1 - ^ - Li/l - 



sty \ s / \ t / V s / V t 

Li 2 (l - aP 2 ) + Li 2 (l - aQ 2 ) - Li 2 (l - as) - Li 2 (l - at) 



1 . 2 / s\ 7T 



r 2 



+ - log" - - + — + U = 



2 ° V tj 2 
where 

M := , og (m i og r i - ^ - io g * u- p -\- * ( p -\ * u - p * 



st J \ st J \ s J \ s J \ t J \ t J 

+ log(aP 2 ) log(l - aP 2 ) + log(aQ 2 ) log(l - aQ 2 ) - log(as) log(l - as) - log(at) log(l - at) . 
Our identity (j5.18|) is proved if T> = 0, where 

P := W + ilo^(f) + I log =(-|) + f (C.17) 

We have checked that the right-hand side of (jC17|) indeed vanishes whenever Mantel's 
identity (jC.lOJ) is satisfied. 
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